>  I 


«r 


DISTRIBUTION  STATEMENT  A 
Approved  ior  Public  Release 
Distribution  Unlimited 


TWO-DIMENSIONAL  HEAT  TRANSFER  IN 
RADIATING  STAINLESS-STEEL-CLAD  COPPER  FINS 


by  Norbert  0.  Stockman  and  Edward  C.  Bittner 
Lewis  Research  Center 


Cleveland,  Ohio 


20060516238 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION  -  WASHINGTON,  D.  C.  -  DECEMBER  1965 


NASA  TN  D-3102 


TWO-DIMENSIONAL  HEAT  TRANSFER  IN  RADIATING 

STAINLESS-STEEL- CLAD  COPPER  FINS 

By  Norbert  O.  Stockman  and  Edward  C.  Bittner 

Lewis  Research  Center 
Cleveland,  Ohio 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


For  sale  by  the  Clearinghouse  for  Federal  Scientific  and  Technical  Information 
Springfield,  Virginia  22151  -  Price  $1.00 


TWO-DIMENSIONAL  HEAT  TRANSFER  IN  RADIATING  STAINLESS-STEEL-CLAD  COPPER  FINS 
by  Norbert  0.  Stockman  and  Edward  C.  Bittner 
Lewis  Research  Center 


SUMMARY 


r* 


An  analysis  of  the  two-dimensional  heat  transfer  in  a  radiating  composite 
fin  is  presented  together  with  details  of  the  numerical  method  of  solution.  ^ 
Several  stainless-steel-clad  copper  fins* typical  of  finned-tube  space  radiators 
are  analyzed  over  a  range  (0  to  l)  of  stainless -steel  to  copper  thickness  ratio. 
Items  discussed  are  temperature  profiles  along  the  fin  length,  temperature  drop 
across  fin  thickness,  effective  conductivity,  fin  efficiency,  heat  rejection 
rate  per  unit  weight,  and  ratio  of  effective  conductivity  to  effective  density. 
It  was  found  for  the  particular  geometries  analyzed  that  the  two-dimensional 
temperature  distribution  in  the  fin  cross  section  is  essentially  one¬ 
dimensional.  The  use  of  the  thickness -averaged  conductivity  in  one-dimensional 
calculations  yields  essentially  the  same  results  as  the  two-dimensional  calcu¬ 
lations  .  * 
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INTRODUCTION 

. . . . — - .  . / 

The  radiator  of  an  electric -power  generation  syspem  for  qpqce  applications 
may  account  :pr  a  substantial  portion  of  the  powerplaht  total 'yfeight  (refs.  1 
and  2).  An  important  factor  in  the  minimum-weight  design  of  these  radiators  is 
the  selection  of  suitable  materials  (ref.  3).  The  large  number  of  basic  re¬ 
quirements  (see  ref.  3  for  discussion)  that  must  be  satisfied  makes  the  selec¬ 
tion  difficult.  Often  all  the  desirable  qualities  cannot  be  found  in  any  one 
currently  available  material,  and  composite  materials  are  therefore  considered. 
In  particular,  the  fin  material  for  a  space  powerplant  radiator  should  have  a 
high  value  for  the  ratio  of  thermal  conductivity  to  density  (ref.  3),  adequate 
strength,  and  a  low  sublimation  rate  in  vacuum  at  high  temperature.  Stainless 
steel  has  adequate  strength  and  satisfactory  sublimation  rate  but  low  thermal 
conductivity,  whereas  copper  has  very  high  conductivity  but  questionable 
strength  and  sublimation  rate  at  temperatures  above  about  800°  K.  These  facts 
lead  to  a  consideration  of  the  composite  material,  stainless-steel-clad  copper, 
for  radiator  fins  (e.g.,  in  the  SNAP-50/SFUR  design  studies,  refs.  4  and  5). 
Stainless-steel-clad  copper  has  already  been  used  in  other  applications  (e.g., 
cooking  utensils  and  heat  exchangers)  for  several  years  and  is  commercially 
available.  If  a  composite  material  is  to  be  used  for  radiator  fins,  its  effec¬ 
tive  thermal  conductivity  must  be  known  for  use  in  one- dimensional  calculations 
for  radiator  design  and  comparison  purposes.  Suppliers  of  stainless -steel-clad 
copper  material  usually  quote  two  values  of  effective  conductivity:  one  value 


for  one -dimensional  conduction  in  the  direction  parallel  to  the  plane  of  the 
interface,  and  another  value  for  one -dimensional  conduction  normal  to  the  plane 
of  the  interface.  For  any  significant  thickness  of  clad  these  two  values  are 
quite  different  (e.g.,  for  a  stainless-steel  to  copper  thickness  ratio  of  only 
0.05  these  values  are,  respectively,  0.955  and  0.525  that  of  copper  alone).  In 
a  radiator  fin  the  conduction  heat  transfer  is  two-dimensional,  and  the  value 
of  effective  conductivity  will  lie  somewhere  between  the  two  one-dimensional 
values . 

Specifically,  the  heat  radiated  by  a  fin  is  initially  supplied  across  the 
fin  thickness  at  the  tube  edge  and  is  conducted  parallel  to  the  fin  surface. 

The  heat  must  subsequently  reach  the  surface  of  the  fin  where  it  is  radiated  to 
space.  There  is  clearly  a  two-dimensional  effect  involved,  and  the  magnitude 
of  this  effect  will  increase  as  the  stainless-steel  to  copper  thickness  ratio 
increases  since  the  stainless  steel  has  a  much  lower  conductivity  than  copper. 
On  the  other  hand,  a  radiator  fin  is  usually  very  long  relative  to  its  thick¬ 
ness,  which  will  tend  to  decrease  the  two-dimensional  effect.  Therefore,  it  is 
not  a  priori  evident  how  great  the  two-dimensional  effect  will  be,  nor  is  it 
evident  what  the  effective  conductivity  of  a  radiating  fin  will  be. 

Reference  6  attacks  this  type  of  problem  in  a  quasi- two-dimensional  way  by 
assuming  a  one -dimensional  heat  flow  in  the  fin  core  parallel  to  the  radiating 
surface  and  one-dimensional  heat  flow  in  the  clad  or  coating  normal  to  the  ra¬ 
diating  surface.  The  primary  application  of  the  analysis  of  reference  6  is  for 
coated  fins  having  a  relatively  thin,  very  low  conductivity,  nonmetallic  coat¬ 
ing.  The  metallic-clad  and  the  nonmetallic-coated  fins  present  essentially  the 
same  problem.  However,  it  is  not  obvious  that  the  method  of  reference  6  will 
be  adequate  for  metallic-clad  fins  in  all  cases  because  a  metallic  clad  is 
thicker  and  has  a  higher  value  of  thermal  conductivity  than  a  nonmetallic  coat¬ 
ing. 


This  report  presents  a  brief  analysis  of  the  two-dimensional  heat  transfer 
in  a  composite  fin  with  the  conventional  fin-tube  radiator  boundary  conditions, 
that  is ,  radiation  from  both  surfaces  and  prescribed  constant  temperature  at 
the  fin  base.  An  exact  effective  conductivity  is  obtained  from  the  two- 
dimensional  results  and  is  compared  with  the  conductivity  for  one-dimensional 
conduction  in  the  direction  parallel  to  the  fin  surface. 

Several  cases  are  analyzed  that  cover  a  range  of  the  conductance  parameter 
( the  conventional  identifying  parameter  for  radiating  fins ) .  The  range  in  this 
parameter  is  obtained  by  varying  the  length  and  thickness  of  the  fin  at  two 
different  temperature  levels.  For  each  case  several  values  of  stainless-steel 
to  copper  thickness  ratio  are  prescribed.  The  results  to  be  presented  are 
discussed  under  two  main  headings:  first,  the  degree  of  two-dimensionality  of 
the  several  cases,  and  second,  the  effect  of  varying  the  conductance  parameter 
and  the  stainless -steel  to  copper  thickness  ratio  on  fin  temperature  profiles, 
fin  efficiency,  heat  rejected  per  unit  weight,  and  effective  conductivity  to 
effective  density  ratio. 
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(b)  As  considered  in  analysis  (thickness  exaggerated). 
Figure  1.  -  Cross  section  of  stainless-steel-clad  copper  fin. 


ANALYSIS 

The  configuration  to  he  analyzed 
is  shown  in  figure  1.  This  is  a 
stainless -steel-clad  copper  fin  such 
as  might  he  used  in  a  central  fin-tuhe 
radiator  (fig.  1(a)).  Only  one  quar¬ 
ter  of  the  fin  is  analyzed  (fig.  1(b)) 
because  symmetry  is  assumed  about 
x  =  L  and  y  =  0.  There  is  no  heat 
flow  across  these  boundaries.  The 
temperature  at  x  =  0  is  assumed  con¬ 
stant  and  equal  to  the  tube  outer- 
surface  (base)  temperature  in  a  fin- 
tube  configuration. 

The  analysis  is  based  on  steady- 
state  two-dimensional  conduction  in 
the  fin  and  the  clad  with  perfect  ther¬ 


mal  contact  at  the  interface  (y  =  tcu). 
The  surface  of  the  clad  (y  =  t)  is 
assumed  to  be  radiating  hemispherically 
at  constant  emittance  to  space,  but  radiant  interchange  with  the  tubes  and  in¬ 
cident  radiation  from  the  surroundings  are  neglected  for  simplicity.  Also,  the 
effective  sink  temperature  of  space  is  assumed  to  be  zero.  It  is  assumed  that 
the  thermal  conductivity  along  the  fin  does  not  vary  with  temperature.  For  the 
range  of  variables  considered  herein  the  maximum  error  in  assuming  constant  con¬ 
ductivity  is  estimated  as  less  than  0.6  percent  by  the  method  of  reference  7. 


Basic  Equations 

The  equation  governing  the  two-dimensional  conduction  in  the  copper  is 

S2Tcu(x,y)  a2Tcu(x,y) 
dx2  dy2 

and  in  the  stainless  steel  clad  is 

a2Tss(x,y)  ^  d2Tss(x,y) 

Sx2  ay2  “  1 

where  Tcu(x,y)  is  the  local  temperature  in  the  fin  and  Tss(x,y)  is  the  local 
temperature  in  the  clad.  (All  symbols  are  defined  in  appendix  A.) 

At  the  boundary  y  =  tcu  between  the  fin  and  clad,  perfect  thermal  con¬ 
tact  is  assumed: 


®cu(x^cu)  ~  ^ss(x,tcu) 


(3) 
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Also  at  this  boundary  there  must  be  continuity  of  heat  flow;  hence. 


In  order  to  facilitate  comparison  of  the  solutions  for  different  cases, 
the  following  dimensionless  parameters  are  formed: 

(1)  Clad-to-core  thickness  ratio,  tBS/tcu 

(2)  Aspect  ratio  of  copper  fin,  tcu/L 

(3)  Conductivity  ratio,  kss/kcu 

(4)  Temperature  ratio,  @(X, Y)  =  T(x,y)/T0 

(5)  Distance  ratio  for  any  distance,  X  =  x/L  and  Y  =  y/L 

(6)  Conductance  parameter,  X  =  eaT^L^/kcut 

This  last  ratio  does  not  occur  naturally  in  the  two-dimensional  analysis  of 
clad  fins,  but  it  is  used  because  of  its  wide  occurrence  in  fin-tube  radiator 
literature . 

Putting  the  dimensionless  parameters  into  equations  (l)  and  (2)  yields 
the  following: 
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In  the  fin: 


In  the  clad: 


At  Y  = 


-cu 


and 


At  Y  =  i=hE|l+  "S3 


L  L 


cu 


At  X  =  0: 


At  X  =  1.0: 


a2ecu(x,Y)  a2ecu(x,Y) 
ax2  sy2 

(6) 

a2ess(x,Y)  |  a2ess(x,Y)  o 

Sx2  bY2 

(7) 

then  become  the  following: 

ecu(x>  -p)  “  ess(x, 

(8) 

aecu(x,Y)  kss  Sess(x,Y) 

5y  kcu  8i 

(9) 

/ 

aess<x’Y>  ,  *cu  kOU  .4  /-,  t\ 

5y -  -  “A  L  k5S  ess^Lj 

(10) 

ess(x,Y)  =  ecu(x,Y)  =  1.0 

(11) 

aess(x^)  „ 

ax  ax 

(12) 

At  Y  =  0 : 


£>0cu(X;Y) 

bY 


=  0 


(13) 


It  is  thus  seen  that  the  independent  parameters  are  teu/L,  tss/tcu,  kss/kcu> 
and  A,  and  their  specification  determines  the  temperature  distribution  through¬ 
out  the  fin.  Results  are  given  in  dimensionless  form  in  terms  of  these  dimen¬ 
sionless  variables;  however,  the  calculations  were  made  on  the  dimensional  eq.ua- 
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tions  (l)  and  (2).  These  equations  subject  to  the  dimensional  boundary  condi¬ 
tions  given  previously  were  solved  numerically  on  an  IBM  7094  computer  using  a 
finite -difference,  line-by-line,  overrelaxation  method.  Details  of  the  numer¬ 
ical  method  of  solution  are  given  in  appendix  B. 


Fin  Efficiency 


A  common  measure  of  fin  performance  is  the  radiating  efficiency,  which  is 
defined  as  the  ratio  of  the  actual  heat  radiated  by  the  fin  based  on  the  calcu¬ 
lated  temperature  distribution  along  the  outer  surface  to  the  heat  radiated  by 
an  isothermal  fin  at  the  base  temperature  T0:  that  is. 


(14) 


This  parameter  will  be  used  to  evaluate  the  effect  of  various  thicknesses  of 
clad. 


Effective  Thermal  Conductivity 

The  effective  conductivity  of  a  clad  fin  is  a  function  of  the  conductiv¬ 
ities  and  the  relative  dimensions  of  the  copper  core  and  the  stainless  steel 
clad.  It  is  also  a  function  of  the  path  the  flow  of  heat  must  take  to  reach 
the  radiating  surface  of  the  clad  fin.  These  considerations  are  reflected  in 
an  effective  conductivity  k  of  clad  fin  which,  in  this  report,  is  that  value 
of  conductivity  that  satisfies 


Aeff  = 


eaT^L2 

kt 


(15) 


where  Aeff  is  that  value  of  A  which,  when  used  in  one -dimensional  radiating 
fin  calculations  (such  as  those  of  refs.  8  and  9  or  the  constant  property  cases 
of  ref.  7),  will  result  in  the  value  of  radiating  efficiency  t|  obtained 
herein  by  two-dimensional  calculations. 


BESULTS  AND  DISCUSSION 

Several  fin  geometries  representative  of  300-kilowatt  and  1-megawatt 
electrical-power-output  Rahkine-cycle  space  powerplants  were  analyzed.  Two 
temperature  levels  (811°  and  950°  K)  were  chosen,  and  the  fin  length  and  fin 
thickness  were  varied  to  obtain  a  range  of  practical  values  of  the  identifying 
parameter  A.  These  input  parameters  are  summarized  in  table  I.  Each  basic 
configuration  of  table  I  was  analyzed  over  a  range  of  thickness  ratios  between 
0  and  1.  The  variation  in  thickness  ratio  for  a  given  configuration  was  ob- 
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Temperature  ratio,  T /T( 


TABLE  I.  -  IHPUT  PARAMETERS8, 


[Emittance,  €,  0.9;  thermal  conductivity  of  copper,  kcu, 
3.65  W/(cm)(°K);  thermal  conductivity  of  stainless 
steel,  kss,  0.1818  W/(cm)(°K);  conductivity  ratio, 

kss/kcu>  °-05-] 


Case 

Base 

temper¬ 

ature, 

To> 

°K 

Fin 

length, 

L, 

cm 

Copper 

fin 

thickness, 

w 

cm 

Conductance 
parameter , 

A 

Aspect 
ratio  of 
copper 
fin, 

t  /l 

cu' 

A 

811 

2.0 

0.05 

0.1218 

0.0250 

B 

5.1 

.20 

.1904 

.0400 

C 

4.1 

.05 

.48  74 

.0125 

D 

10.2 

.20 

.7616 

.0200 

E 

5.6 

.05 

1.0000 

.0089 

F 

15.3 

.20 

1.7136 

.0131 

G 

950 

1.8 

0.05 

0.1499 

0.0286 

H 

4.3 

.20 

.2210 

.0471 

I 

3.6 

.05 

.5996 

.0143 

J 

8.6 

.20 

.8841 

.0236 

K 

12.9 

.20 

1.9892 

.0155 

aTo  convert  from  cm  to  in.,  multiply  by  0.3937;  °K  to  °F, 

°F  =  9/5  °K  -  460;  °K  to  °R,  multiply  by  9/5;  W/(cm)(°K)  to 
Btu/(hr )(ft)(^R),  multiply  by  57.8. 


tained  by  decreasing  tcu,  the 
thickness  of  copper,  and  in¬ 
creasing  tss,  the  thickness  of 
stainless  steel  clad,  while  main¬ 
taining  a  constant  total  thick¬ 
ness  t.  The  conductivity  ratio 
kss/kcu  should  be  a  function  of 
fin-base  temperature;  however, 
to  reduce  the  number  of  parame¬ 
ters,  the  same  value  was  used  for 
both  temperature  levels. 

The  results  to  be  presented 
are  discussed  under  two  main 
headings:  first,  the  degree  of 
two  dimensionality  present  in 
the  cases  analyzed,  and  second, 
the  effect  on  certain  fin  parame¬ 
ters  of  varying  the  stainless - 
steel  to  copper  thickness  ratio. 


Figure  2.  -  Surface  temperature  distribution  of  cases  A  and  K  for  thickness 
ratios  from  0  to  1. 0, 


Degree  of  Two-Dimensionality 


The  items  to  be  discussed 
in  this  section  are  temperature 
profiles  along  the  fin  length, 
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temperature  drop  across  the  fin 
thickness ,  effective  conductiv¬ 
ity,  and  comparison  of  two- 
dimensional  and  one-dimensional 
results . 

Temperature  drop  across  fin.  - 
Of  the  several  cases  analyzed 
(table  I),  case  A  has  the  smallest 
temperature  gradient  along  the  fin 
length  and  case  K  has  the  largest. 
The  temperature  profiles  of  these 
two  extreme  cases  at  three  values 
of  thickness  ratio  tss/tcu  are 
shown  in  figure  2.  The  two- 
dimensional  effect,  that  is,  the 
difference  between  the  tempera¬ 
ture  profile  along  the  radiating 
surface  of  the  fin  (Y  =  t/L)  and 
that  along  the  adiabatic  surface 
(Y  =  0),  is  not  large  enough  to 
be  seen  on  the  scale  of  figure  2 . 


Figure  3.  -  Temperature  ratio  drop  across  fin  thickness  against  dis¬ 
tance  ratio  along  fin. 

figure  3  in  the  form  of  the  temperature  ratio  drop  across  the  fin  from  *Y  =  0 
to  Y  =  t/L  for  three  values  of  tss/tcu.  The  cases  shown  are  case  K,  because 
it  has  the  steepest  X-direction  temperature  gradient  (fig.  2),  and  case  H,  be¬ 
cause  it  has  the  largest  Y-direction  temperature  gradient.  The  variation  in 
^ss/^cu  was  obtained  by  reducing  tQU  and  increasing  t5S  to  maintain  a  con¬ 
stant  total  thickness  t.  The  clad  fins  have  a  greater  two-dimensional  effect 
than  the  all-copper  fins  (tss/tcu  =  0)  as  shown  in  the  figure.  It  should  be 
noted,  however ,  that  the  scale  of  the  ordinate  in  figure  3  is  greatly  expanded 
and  the  maximum  difference  in  temperature  ratio  between  the  fin  centerline  and 
the  fin  surface  is  only  about  1.15x10" 3  unit. 


The  two-dimensional  effect 
is  shown  more  specifically  on 


Effective  conductivity.  -  As  was  stated  in  the  INTRODUCTION  the  value  of 
effective  conductivity  of  a  radiating  clad  fin  should  lie  somewhere  between  the 
one-dimensional  value  for  conduction  in  the  direction  parallel  to  the  fin  sur¬ 
face  and  that  for  conduction  normal  to  the  surface.  In  view  of  the  small  two- 
dimensional  effect  noted  in  the  previous  section,  it  is  to  be  expected  that  the 
effective  conductivity  of  the  radiating  fin  would  be  close  to  that  for  one¬ 
dimensional  conduction  in  the  parallel  or  x-direction.  To  check  this  expecta¬ 
tion,  the  exact  effective  conductivity  for  the  two-dimensional  radiating  fin 
will  be  investigated  and  will  then  be  compared  with  the  one-dimensional 
x-direction  conductivity. 


An  exact  value  of  effective  conductivity  can  be  obtained  from  equation  (15). 
Solving  equation  (15)  for  the  exact  effective  conductivity  k  and  forming  the 
ratio  k/kcu  yield 
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3  -9 

|JC 


X  =  £°t°l2  (16' 

Au  ^AuAff 

The  grouping  eaTolP/tkcll  of  equation 
(16)  is  equal  to  the  conductance  parame¬ 
ter  A.  Substituting  A  into  equa¬ 
tion  (16)  yields 


^5^.  A  value  of  Aepp  was  obtained  from 

one -dimensional  calculations  for  each 
.  5 _ 1 _ _ 1 _ _ _ i  value  of  t]  obtained  from  the  two- 

0  *(flad-to-cor"e* thickness" ratio,  t,sft'cf.  l#  dimensional  calculations  reported  herein, 

and  the  effective  conductivity  ratio  was 
Figure  4  -  Exact  effective  conductivity  ratio  against  thick-  obtained  from  equation  (17).  On  fig- 

ure  4  are  shown  curves  for  the  exact  ra- 
_  diating-fin  effective  conductivity  ratio 

k/kcu  as  a  function  of  thickness  ratio  tss/tcu  for  cases  E  and  H.  The 
remaining  cases  listed  in  table  I  fall  between  these  two  curves  and  are  omitted 
for  clarity.  The  close  agreement  among  the  several  cases  indicates  that  the 
effective  conductivity  ratio  k/kcu  is  dependent  primarily  on  thickness  ratio 
AsAcu  for  a  given  conductivity  ratio  kss/kcu.  The  exact  value  of  effective 
conductivity  will  be  compared  with  the  effective  conductivity  based  on  one¬ 
dimensional  conduction  in  the  x-direction,  which  is  the  thickness-averaged 
value  given  by 


Askss  +  Aukcu 

As  +  Au 


Writing  equation  (18)  in  terms  of  dimensionless  ratios  yields 


As  As 
Au  kcu 


Thus,  the  ratio  kx/kcu  is  a  function  only  of  the  thickness  ratio  iss/tcu 
and  the  conductivity  ratio  kss/kcu.  The  relation  between  K^/keu  and  tss/tcu 
is^ shown  on  figure  5(a)  for  a  conductivity  ratio  kss/kcu  of  0.05.  Also 
plotted  in  figure  5(a)  is  the  exact  effective  conductivity  ratio  of  the  case 
that  deviates  most  from  the  thickness-averaged  conductivity  (case  H) .  The  close 
agreement  between  kx/kcu  and  k/kcu  is  evident. 


As  was  stated  earlier,  only  one  value  of  kss/kcu  was  used  even  though 
two  temperature  levels  were  considered.  To  determine  whether  the  change  in 
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(a)  Conductivity  ratio,  0.05.  (b)  Conductivity  ratio,  0.07. 


Figure  5.  -  Comparison  of  thickness-averaged  conductivity  with  exact  two-dimensional  effective  conductivity. 


kss/kcu  over  the  temperature  range  of  interest  significantly  affected  the 
agreement  of  k/kcU  and  kx/kCu>  case  H  was  analyzed  with  a  conductivity  ratio 
kss/kcu  0*°7  (corresponding  to  a  temperature  of  911°  K) .  The  results  are 
shown  in  figure  5(b),  and  the  agreement  is  seen  to  be  not  significantly  differ¬ 
ent  from  that  of  figure  5(a). 

This  close  agreement  of  kx/kcu  with  k/kcu  means  that  the  thickness- 
averaged  effective  conductivity  can  safely  be  used  in  radiator  heat-transfer 
calculations  for  geometries  comparable  to  those  treated  herein. 

Comparison  of  one -dimensional  and  two-dimensional  results.  -  As  a  further 
check  of  the  closeness  between  the  results  of  one -dimensional  calculations  with 
the  thickness-averaged  conductivity  and  the  two-dimensional  solutions, 

several  cases  were  analyzed  by  the  method  of  reference  7  but  with  constant 
properties.  (The  one -dimensional  methods  of  refs.  8  or  9  could  also  be  used.) 
The  cases  chosen  from  table  I  were  those  having  extreme  values  of  A,  cases  A 
and  K,  and  those  having  extreme  values  of  t/L,  cases  E  and  H  (case  H  has  the 
greatest  two-dimensional  effect).  In  all  cases  the  temperature  profiles  and 
fin  efficiencies  agreed  with  the  two-dimensional  results  to  within  three  places 
with  small  variations  in  the  fourth  place.  Thus,  for  the  range  of  A,  t/L,  and 
tss/tcu  treated  herein,  the  use  of  thickness -averaged  conductivity  in  one¬ 
dimensional  calculations  will  give  essentially  the  same  results  as  the  more 
exact  two-dimensional  calculations. 


Effect  of  Varying  Thickness  Ratio 

The  stainless-steel  to  copper  thickness  ratio  used  will  in  general  depend 
on  the  particular  application.  All  the  cases  of  table  I  were  analyzed  with 
thickness  ratios  tss/tcu  of  0,  0.1,  0.3,  0.5,  and  1.0.  The  effect  of  this 
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Heat  rejection  per  unit  weight  ratio,  (Q/W)/(Q/W) 


Figure  6.  -  Fin  efficiency  against  thickness  ratio  for  cases  with  constant 
total  fin  thickness. 


Figure  7.  -  Heat  rejection  per  unit  weight  ratio  against  thickness  ratio. 


variation  in  thicknsss  ratio  on 
several  parameters  used  in  evalu¬ 
ating  fins  or  fin  materials  is 
presented  in  this  section.  Most 
of  the  results  presented  were 
obtained  from  the  same  two- 
dimensional  calculations  as  the 
temperature  profiles  presented 
previously.  Essentially  the 
same  results  could  also  have 
been  obtained  from  one¬ 
dimensional  calculations  because 
of  the  close  agreement  noted  in 
the  preceding  section. 

Fin  efficiency.  -  The  effect 
on  fin  efficiency  T)  of  varying 
tgs/tcu  is  shown  on  figure  6. 
This  figure  indicates  a  signifi¬ 
cant  drop  in  heat  transfer  as 
the  high  conductivity  copper  is 
replaced  by  the  lower  conductiv¬ 
ity  steel.  Thus,  making  the 
simplest  assumption  and  treating 
the  clad  as  if  it  were  copper 
would  introduce  a  large  heat  re¬ 
jection  error  in  radiator  calcu¬ 
lations  (except  for  very  small 
tss/tcu)*  For  example,  for  case  K 
(A  =  1.99)  with  tss/tcu  =  0.5, 
the  fin  efficiency  is  16  percent 
lower  than  that  for  the  all 
copper  fin  (tss/tcu  =  0).  How¬ 
ever,  as  was  indicated  earlier, 
the  clad  fin  can  be  treated 
accurately  in  one -dimensional 
radiator  calculations  with  the 
use  of  the  thickness -averaged 
( effective )  conductivity . 

Heat  rejection  rate  per  unit 
weight.  -  An  important  parameter 
used  for  comparative  evaluation 
of  a  radiator  component  is  the 
heat  rejection  rate  per  unit 
weight  Q/W.  Therefore,  to  study 
the  effect  of  varying  the  propor¬ 
tion  of  stainless  steel  to  copper 
in  a  clad  fin  the  ratio 
(q/w)/(q/w)  is  formed,  where 
(Q/w)0  is  the  Q/w  for  an  all 
copper  fin  (tss/tcu  =  0).  This 
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Figure  8.  -  Ratio  of  effective  conductivity  to  effective  density 
against  thickness  ratio. 


stainless -steel-clad  copper  on  this  ba: 
is  the  average  density  of  the  clad  fin 


ratio  is  plotted  in  figure  7  as  a 
function  of  tss/tcu.  For  most  practi 
cal  radiators,  tss/tcu  will  be  less 
than  0.5,  and  consequently  the  reduc¬ 
tion  in  Q/W  will  vary  between  3 
and  14  percent  for  values  of  A 
ranging  between  0.12  and  2.0,  respec¬ 
tively.  However,  this  is  the  effect 
on  the  fin  only.  The  net  effect  of 
using  clad  fins  on  the  total  radiator 
Q/w  will  depend  on  the  relative  con¬ 
tributions  of  the  fins  to  total  radia¬ 
tor  weight  and  in  many  cases  will  be 
less  than  the  effect  shown  in  figure  7 

Conductivity  to  density  ratio.  -  A 
parameter  based  on  material  properties 
that  is  used  to  compare  the  conducting 
capabilities  of  materials  on  a  weight 

basis  is  k/p_(ref.  3).  To  compare  . 

,  the  ratio  kx/p"  is  formed,  where  p 
Ld  is  given  by 


P 


^ssPss  +  ^cu^cu 
tss  +  tcu 


(20) 


and  kx  is  the  thickness -averaged  conductivity  defined  previously  in  equa¬ 
tion  ( 18 ) . 


TABLE  II.  -  C ONDU CTTV ITY  TO  DENSITY  RATIOa  FOR 
SEVERAL  MATERIALS  AT  811°  K 


Material 

Conductivity  to 
density  ratio, 

(w)(m2,/(kg)(°K) 

Pyrolytic  graphite 

0.132 

Beryllium 

.0583 

Commercial  graphite 

.0550 

Copper,  tss/tcu  =  0 

.0428 

Stainless-steel-clad  copper,  tsfl/tnn  =0.5 

.0308 

Stainless-steel-clad  copper,  tss/tcu  =  1.0 

.0231 

Molybdenum 

.0132 

Columbium 

.0077 

Tantalum 

.00385 

Titanium 

.00385 

Stainless  steel 

.00286 

Vanadium 

.00231 

aAll  data  except  stainless-steel-clad  copper  obtained  from 
fig.  9  of  ref.  3. 
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The  ratio  E^/p  is  plotted  against  tss/tcu  in  figure  8.  In  table  II 

these  values  of  kx/ p  are  compared  with  the  k/p  values  obtained  for  a  variety 
of  materials  from  reference  3.  It  can  be  seen  that  for  thickness  ratios  up  to 
0.5,  stainless-steel-clad  copper  has  a  higher  k/p  value  than  most  of  the 
materials  listed,  and  is  not  too  far  removed  from  the  best  materials. 


SUMMARY  OP  RESULTS 

For  the  particular  configurations  analyzed  the  following  results  have  been 
obtained : 

1.  The  temperature  in  the  clad  fin  shows  only  a  slight  two-dimensional 
effect. 


2.  The  effective  conductivity  of  a  radiating  clad  fin  is  very  nearly  equal 
to  the  thickness-averaged  conductivity  (one-dimensional)  of  the  clad  fin. 

3.  One -dimensional  calculations  that  use  the  thickness-averaged  conductiv¬ 
ity  yield  essentially  the  same  results  for  temperature  profile  and  fin  effec¬ 
tiveness  as  the  two-dimensional  calculations. 

4.  Increasing  the  proportion  of  stainless  steel  in  a  clad  fin  of  given 
total  thickness  decreases  the  radiating  effectiveness  of  the  fin.  For  a 
stainless-steel  to  copper  thickness  ratio  of  0.5,  the  reduction  is  16  percent. 

5.  The  ratio  of  effective  conductivity  to  effective  density  of  a  stainless- 
steel-clad  copper  fin  is  relatively  high  and  remains  competitive  with  the  other 
materials  of  interest. 


Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  August  17,  1965. 
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APPENDIX  A 


s* 

a 


b 

<4 

c 

— > 
d 

I 

J 

k 

k 

L 

m 

NDX 

NDYC 

NDYF 

NX 

NY 

n 

Q 

^ideal 


SYMBOLS 

matrix  of  system  of  finite-difference  equations 
diagonal  elements  of  matrix  Si 

submatrices  of  order  NY  of  si  which  lie  along  the  diagonal 
off-diagonal  elements  of  matrix  2 
off-diagonal  submatrix  of  order  NY  of  si 
elements  of  matrix  # 

vector  whose  components  are  right-hand  sides  of  finite -difference 
equation 

position  number  of  horizontal  grid  line 
position  number  of  vertical  grid  line 
thermal  conductivity,  W/(cm)(°K) 

effective  conductivity  of  radiating  clad  fin,  W/(cm)(°K) 

effective  conductivity  of  clad  fin  transferring  heat  by  conduction  in 
x-direction  only,  W/(cm)(°K) 

fin  length,  cm 

number  of  current  iteration 

number  of  vertical  grid  strips 

number  of  horizontal  grid  strips  in  clad 

number  of  horizontal  grid  strips  in  fin 

total  number  of  vertical  grid  lines 

total  number  of  horizontal  grid  lines 

distance  normal  to  surface,  cm 

heat  rejection  rate,  W/(cm)(sec) 

heat  rejection  rate  of  an  isothermal  fin  at  temperature  T  ,  . 
w/(cm)(sec) 
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s  distance  along  a  surface,  cm 

T  temperature,  °K 

T  vector  of  temperatures 

— > 

T  vector  of  temperatures,  intermediate  in  the  computation 

t  thickness,  cm 

¥  fin  weight,  kg 

X  x/L 

x  distance  along  length  of  fin,  cm 

Ax  grid  spacing  in  x-direction,  cm 

Y  y/L 

y  distance  along  thickness  of  fin,  cm 

£y  grid  spacing  in  y-direction,  cm 

z  dummy  variable  vector 

Zj:  components  of  z 

a-j_  coupling  parameters  (eg..  (B15))  in  difference  equations 
e  emittance 

r)  radiating  efficiency  of  fin 

6  temperature  ratio,  T/T0 

3  2/ 

A  conductance  parameter,  €0ToL  /kcut 

Aeff  A  of  clad  fin,  eaT^L^/kt 

p  vector  of  ratio  of  elements  of  the  matrix  £8 

p^  components  of  p 

p  density,  kg/m® 

p  thickness -averaged  density,  kg/m® 

0  Stefan-Boltzmann  constant,  5.67X10"®  ¥/(m®)(°K^) 

co  block-successive  overrelaxation  parameter 
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CD]-,  optimum  value  of  cd 

Subscripts : 
cu  copper  fin 

J  variable  horizontal  position  along  fin 

\ 

NX  last  horizontal  position  along  fin 

NY  last  vertical  position  on  fin 

o  value  at  x  =  0  or  tss/tcu  =  0 

ss  stainless  steel  clad 

Superscripts : 

(m)  m"kh  iterate 

(o)  initial  approximation 
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APPENDIX  B 


NUMERICAL  APPROXIMATION  OF  TEMPERATURE  DISTRIBUTION 
Equations  (l)  and  (2)  may  be  replaced  by  the  single  equation 


which  holds  in  the  entire  fin  as  a  partial  differential  equation  with  discon¬ 
tinuous  coefficients.  Such  problems  have  been  treated  recently  (e.g.,  ref.  10). 
To  solve  equation  (Bl)  subject  to  the  boundary  conditions,  a  rectangular  grid 
is  superimposed  on  the  composite  rectangle  above  the  x-axis  of  figure  l(b)  (p.3) 
in  such  a  way  that  each  boundary  (including  the  interface)  is  a  grid  line.  Ax 
is  a  constant  L/NDX  where  NDX  is  the  number  of  vertical  strips  laid  off  in 
the  rectangle  by  the  discretization,  and  Ay  is  one  constant  in  the  fin 
Aycu  =  tcu/NDYF  an<3-  another  in  the  clad  Ayss  =  tss/NDYC.  The  grid  points  are 
labeled  with  ordered  pairs  (l,J)  of  positive  integers  in  the  manner  to  be  de¬ 
scribed. 

Let  NY  =  NDYF  +  NDYC  +  1  be  the  total  number  of  horizontals  in  the  grid. 
Number  these  1,  2  ...  NY  starting  at  the  lowermost  (y=0).  Let  NX  =  NDX  +  1 
be  the  number  of  verticals  in  the  grid.  Number  these  1,  2  ...  NX  from  left 
to  right  starting  at  the  leftmost  (x  =  0).  The  first  element  of  the  ordered 
pair  (l,J)  is  the  position  number  of  the  horizontal  and  the  second  element  is 
that  of  the  vertical;  thus,  (I,J)  is  such  that  1  <  I  <  NY  and  1  <  J  <  NX. 

Next  the  partial  differential  equations  and  boundary  conditions  are  re¬ 
placed  by  a  system  of  algebraic  equations  (some  of  which  are  nonlinear)  ob¬ 
tained  by  integrating  the  partial  differential  equation  over  the  cells  of  the 
"dual  grid"  (ref.  11,  p.  ll),  by  applying  the  divergence  theorem,  and  by  re¬ 
placing  the  normal  derivatives  occurring  in  the  boundary  integrals  by  differ¬ 
ence  quotients,  zero  or  -eaT^,  depending  on  the  position  of  the  cell  in  ques¬ 
tion  (ref.  12,  sec.  6.3). 

This  process  is  illustrated  by  the  following  example  taken  at  the  point 
x  =  Ax,  y  =  t  (i.e.,  (I,J)  =  (NY,2)).  The  cell  in  question  of  the  dual  grid 
(shaded  in  fig.  9)  has  its  boundaries  (except  the  upper  one)  midway  between 
grid  lines  of  the  original  grid.  The  algebraic  equation  at  (NY, 2)  is  obtained 
in  the  following  steps. 

0  “  ff  [is  (*  H)  +  |f  (k  If) dy 

Cell 

0-4  k  ds 

J Around  °n 

cell 
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Figure  9.  -  Detail  of  grid  for  numerical  solution. 


0  - 


^  lr 

“T0  -  T(NY, 2 )~ 

A^ss  +  k 

T(NY  - 

1,2)  -  T(NY,2)" 

Kss 

Ax 

2  ss 

AYss 

Ax 


+  kss 


T(NY,5)  -  T(NY,2 ) 


Ax 


A^ss 


-  ecrT4(NY,2)Ax  ( B2 ) 


This  particular  point  was  chosen  to  illustrate  the  fact  that  sometimes  a  path 
length  in  the  integration  is  only  one-half  a  total  mesh  width  (Ztyss/2)  and  to 
illustrate  the  ease  of  handling  the  two  types  of  boundary  conditions  inherent 
in  this  method  of  discretization  (note  first  and  fourth  terms). 

Repeating  the  integration  on  all  cells  that  surround  the  points  (l,J), 

1  <  I  <  NY  and  2  <  J  <  NX,  yields  a  system  of  NY(NX  -  l)  equations  (of  which 
NX  -  1  are  nonlinear,  viz.,  (NY,J),  2  <  J  <  NX)  in  as  many  unknowns  T(l,J). 

It  is  noted  that  T(l,l)  =  T0  when  !<,!<,  NY. 

A  direct  solution  of  this  system  is  impossible  because  of  the  nonlinearity; 
therefore,  an  iterative  procedure  must  be  used.  Let  T(m)  (I,J)  represent  the 
mth  iterate  for  T(I,J)  where  T(°)  (I,J)  is  some  initial  approximation  to 
T(l,J).  Equations  at  the  points  (NY,J),  2  <  J  <  NX,  are  linearized  in  the  same 
way  that  Newton's  method  for  finding  the  roots  of  equations  is  derived  (see 
ref.  13,  p.  192).  The  term  [T(m)(NY,J)]4  is  replaced  by  the  first  two  terms 
in  its  Taylor  expansion  about  T(m”l  )(NY,J),  that  is  by 

[T(m-1)(NY,J)]4  +  4[T(m"1)(NY,j)]3[T(ra)(NY,j)  -  T(m"1)(NY, J)] 
which  reduces  to 

4[T(m-l)(NY,j)]3T(m)(NY,J)  -  3[T^m_1^(NY, J)]4 

The  resulting  system  of  equations  linear  in  t(“0  can  be  written  in  matrix  nota¬ 
tion  as 

^(m)  =  3  (B3) 

where 

^^[tU-1)]  (E4) 

is  a  square  matrix  of  order  NY(NX  -  l)  and 
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d  =  d(T^m_1))  (B5) 

is  a  column  vector  with  NY(NX  -  l)  components.  ^  is  a  block  tridiagonal 
matrix,  that  is. 


fa 2 

0 

...  0 

0  \ 

/  -**3 

-<e 

...  0 

0  \ 

...  0 

0 

• 

•  •  o 

•  •  o 

0 

•  •  •  • 

•  •  •  • 

•  •  •  “%X-1 

• 

• 

-y 

\  0  0 

0 

...  -Sf 

%x/ 

(B6) 


Here 


Jfj  =  J1^ 111-1 )(  NY,  J)) 


2  <  J  <  NX 


(B7) 


are  square,  tridiagonal  matrices  of  order  NY,  all  of  whose  elements  are  con¬ 
stant  (independent  of  both  J  and  T(ffl"l))  except  for  the  element  in  the  NY, NY 
position  which  depends  on  T(m“l)(NY, J) .  The  remaining  matrix  on  the  diagonal 
of  is 


^NX  =  |  J,(^ra"1)(NY,NX))  (B8) 

and  ^  is  a  diagonal  matrix  (constant)  of  order  NY.  The  iterative  method  may 
be  described  as  a  line-by-line  (vertical)  overrelaxation  scheme,  which  is 
carried  out  in  the  following  way. 


An  initial  approximation  to  the  temperature  is  taken  at  all  points  (e.g., 
t(°)(I,J)  =  T0).  The  temperatures  are  changed  by  sweeping  through  them  from 
left  to  right  an  entire  vertical  at  one  time;  the  sweep  is  repeated  if  the  con¬ 
vergence  criterion  has  not  been  met.  If  attention  is  focused  on  the  verti¬ 
cal,  the  change  is  carried  out  by  the  following  formulas.  Let 


T^(1,J) 


Jm)(NY,J) 


(B9) 


be  the  vector  of  temperatures  along  the  vertical  at  the  m^B  iterate.  It  is 
assumed  that  T^m.|  has  already  been  computed  and  T^.m)  defined  by 


(BIO) 
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Start 


where  0Sj  is  stored  in  computer 
memory  in  partially  factored  form 
since  is  tridiagonal  and 

almost  constant. w  Therefore  the 

— y' 

computation  of  Tj  is  only  a  mat¬ 
ter  of  a  correction  to  the  right- 
hand  side  and  a  back  substitution. 
Then, 


m(m)  =  Z(m-l)  + 
j  J 


CD 


?(m)  _  5(_m-ljj 
(Bll) 


where  cd  is  the  "block  successive 
overrelaxation  parameter"  (ref.  12, 
ch.  4)  that  controls  the  rate  of 
convergence,  and 


0  <  cd  <  2 


(B12) 


To  further  clarify  this,  the 
computation  ( eq.  (BIO))  on  the 
second  vertical  (J  =  2)  is  examined 
in  detail,  and  reference  is  made 
to  the  block  diagram  of  the  com¬ 
puter  program  (fig.  10 ). 


Figure  10.  -  Computer  flow  diagram. 


For  J  =  2, 
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#  =  diag(cj 
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aFY-2 

"bNY-2 
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bNY-l 


0 

0 

0 


-b 


’NY-1 
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(B13) 


(B14) 
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For  convenience,  define 


^ss^-ss 

Ax 

^cu^cu 


a2  Ax 
Axk 


a*  = 


ss 


5  ^ss 
Axk, 


a 


cu 


4  47, 


cu 


j 


(B15) 


Then,  when  equation  (B2)  and  analogous  equations  for  the  other  points  on 
this  vertical  are  considered,  the  elements  of  an&  ^  can  be  written  as 
follows : 


a 


i 


ra2  + 
2  ( cxg 

4  ai  + 

2(a-1 

k.al  + 


a4 

+  a4) 

a2  +  a3  +  a4 
+  a3) 

a3  +  4ecr  AxJtU-^NY,  j)]  3 


i  =  1 

i  -  2  .  .  .  NDYF 
i  >  NDYF  +  1 

i  =  NDYF  +  2  .  .  .  NY  -  1 
i  =  NY 


b . 


l 


fa4 

\a3 


i  =  1  .  .  .  NDYF  -  1 
i  =  NDYF  .  .  .  NY  -  1 


r  (B16) 


a2 

a, 


/2 


=  < 


(aj_  +  ag)/2 


a. 


1 

W2 


i  =  1 

i  =  2  .  .  .  NDYF 
i  =  NDYF  +  1 

i  =  NDYF  +  2  .  .  .  NY  -  1 
i  =  NY 

J 


Therefore,  with  equation  (BIO)  rewritten  for  J  =  2  as 


=9f\r(1m')  +  ^m-D] 

the  computation  for  T^m^  is  carried  out  by  solving  the  following  system  of 
equations  where,  for  brevity,  the  vector 
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(B17) 


replaces  the  vector  Tgm^ : 


alzl  "  b2z2 

+  a2z2  "  b2z3 
“  b2z2  +  a3z3  -  b3z4 


=  C]_[T0  +  T(l,3)] 
=  c2[Td  +  T(2,3)] 
=  c3[Td  +  T(3,3)] 


V  (B18) 


"bWY-2zNY-2  +  aNY-lzNY-l  “  bNY-lzNY 
-  bNY-lzHY-l  +  aNYzNY 


=  c  uy-ilX,  +  T(  NY-1, 3)] 
=  Cuy[T0  +  T(NY,3)] 


Here  the  superscripts  on  T(l,J)  have  been  dropped  to  be  consistent  with  what 
actually  takes  place  in  the  computer  program,  namely,  the  immediate  replacement 

of  by  T^-m^  once  it  is  computed. 


Of  all  the  numbers  aj_,  b^  and  c^,  only  a^y  depends  on  the  fact  that 
computations  are  being  carried  out  on  the  second  vertical  and  at  the  mbb 
iterate.  Thus,  if  some  other  vertical,  other  than  the  last,  were  examined,  the 
only  changes  that  would  occur  in  the  coefficients  in  (B18)  are  as  follows: 

(1)  ajjy  is  replaced  by  its  correct  value. 

(2)  T0  +  T(I,3)  is  replaced  by  T(l,J  -  l)  +  T(l,J  +  l),  I  =  1  .  .  .  NY. 

Therefore,  if  the  previous  system  is  solved  by  the  elimination  of  unknowns,  it 
is  seen  that  as  one  proceeds  from  vertical  to  vertical  and  from  iteration  to 
iteration,  the  same  arithmetic  operations  are  performed  on  the  same  members  with 
the  same  ratios  (pi)  appearing  repeatedly  (at  least  so  far  as  the  left-hand 
sides  are  concerned).  Even  for  the  right-hand  sides  the  same  operations  are 
being  performed  with  the  same  ratios  (pj_)  as  with  the  left-hand  sides,  but  the 
operations  are  performed  on  numbers  that  depend  on  J  and  m.  It  is  therefore 
advantageous  to  compute  and  store  these  pq  ratios  before  any  iterations  are 
performed.  It  is  this  saving  of  the  Pj_  along  with  the  modified  a-j_  that  has 
been  referred  to  previously  by  "triangularizing  of  ^2"  (block  2  of  fig.  10). 
Applying  the  previously  mentioned  operations  and  p.j_  ratios  to  the  iteration- 
and  vertical-dependent  right-hand  sides  (blocks  4  and  5  of  fig.  10)  leads  to  a 
reduced  system  of  equations,  the  last  of  which  involves  only  zjy  (block  6). 

This  is  then  substituted  into  the  preceding  (modified)  equation,  which  involves 
it  with  an<3-  bbe  process  continues  backwards  (in  I)  to  zj_  (blocks  7 

and  8).  The  values  for  the  temperatures  on  this  second  vertical  are  then  re¬ 
placed  by  the  new  values  (block  9): 
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T(I,2)  +  m[Zl  -  T(l,2)] 


1  =  1.  .  .MI 


Numerical  experiments  at  Lewis  indicate  that  there  is  an  optimum  value  cdj-, 
for  cd,  as  in  the  case  of  equations  that  are  all  linear,  where  an  a  priori 
value  can  be  given.  For  example,  for  a  case  where  NDYF  =  4,  NDYC  =  2, 

NDX  =  100,  L  =  0.1333,  tcu  =  tss  =  0.0004167,  kss  =  kcu  =  210,  T0  =  2337, 
e  =  0.9,  and  a  =  0.1713X10-8,  it  is  found  that  cd^  near  1.96  gives  a  relative 
residual  error 


heat  in  -  heat  out 
heat  out 


near  0.0001  in  less  than  300  iterations.  This  is  a  practicable  convergence 
criterion  for  the  problem  at  hand  and  is  so  used  in  the  computer  program 
( block  11  of  fig .  10 ) . 

Here 

k  ds 
on 

end 


heat  in  = 


-L 


Left 


heat  out  = 


Radiating 

boundary 


dT  . 

V"  ds 
on 


Trapezoidal  integration  is  used  to  be  consistent  with  the  finite -difference 
equations.  These  experiments  also  indicate  that  for  co  <  heat  in  increases 
monotonically  to  meet  heat  out,  which  is  decreasing  as  m  increases  if 

t(°)(i,j)  is  chosen  to  be  TQ  everywhere. 

The  program  was  converted  to  double  precision  after  it  was  found  that  the 
large  aspect  ratio  of  the  cells  of  the  grid  led  to  appreciable  loss  of  signifi¬ 
cance.  The  program  takes  less  than  1  minute  to  carry  out  400  iterations  for 
cases  similar  to  those  mentioned  previously. 
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